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ABSTRACT 

Hubble constant Hp and weighted amplitude of matter fluctuations Sg determinations are biased to higher and lower values, 
respectively, in the late universe with respect to early universe values inferred by the Planck collaboration within flat ACDM 
cosmology. If these anomalies are physical, that is, not due to systematics, they naively suggest that Ho decreases and Sg increases 
with effective redshift. Here, subjecting matter density today Qm to a prior, corresponding to a combination of Planck CMB 
and BAO data, we perform a consistency test of the Planck-ACDM cosmology and show that Sg determinations from fos (z) 
constraints increase with effective redshift. Due to the redshift evolution, a ~3o tension in the Sg parameter with Planck at lower 
redshifts remarkably becomes consistent with Planck within lo at high redshifts. This provides corroborating support for an 
Sg discrepancy that is physical in origin. We further confirm that the flat ACDM model is preferred over a theoretically ad hoc 
model with a jump in Sg at a given redshift. In the absence of the CMB+BAO Q, prior, we find that >3o0 tensions with Planck 
in low-redshift data are ameliorated by shifts in the parameters in high-redshift data. Results here and elsewhere suggest that the 
ACDM cosmological parameters are redshift dependent. Fitting parameters that evolve with redshift is a recognizable hallmark 


of model breakdown. 
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1 INTRODUCTION 


Modern cosmology is a pursuit where one theoretically guesses a 
function (or model) and compares it to another a priori unknown 
function, the Hubble parameter H(z), extracted observationally from 
nature. Given the possibilities, the guess is invariably wrong, but 
existing cosmological data can mask disagreement until precision 
improves. Viewed historically, concordance Lambda cold dark 
matter (ACDM) cosmology emerged when a spatially flat, homo- 
geneous, and isotropic universe, filled exclusively with radiation 
and (pressure-less) matter, required a dark energy sector to explain 
observations Riess et al. (1998), Perlmutter et al. (1999). Further 
improvements in data in recent years have unveiled anomalies in the 
Hubble constant Ho := H(z = 0) Aghanim et al. (2020), Riess et al. 
(2022), Freedman (2021), Pesce et al. (2020), Blakeslee et al. (2021), 
Kourkchi et al. (2020), the weighted amplitude of matter fluctuations 
Sg := 6g./Qn,/0.3 Aghanim et al. (2020), Heymans et al. (2013), 
Joudaki et al. (2017), Troxel et al. (2018), Hikage et al. (2019), 
Asgari et al. (2021), Abbott et al. (2022), the lensing parameter 
Alens and/or curvature Q,; Aghanim et al. (2020), Addison et al. 
(2016), Handley (2021), Di Valentino, Melchiorri & Silk (2019), the 
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late-time integrated Sachs—Wolfe (ISW) effect Granett, Neyrinck & 
Szapudi (2008), Kovacs et al. (2017), Kovacs et al. (2019), Kovacs 
et al. (2022), and high-redshift galaxies that seemingly defy ACDM 
expectations Adams et al. (2022), Labbé et al. (2023), Castellano 
et al. (2022), Naidu et al. (2022). See Di Valentino et al. (2021), 
Perivolaropoulos & Skara (2022), Abdalla et al. (2022) for reviews of 
ACDM anomalies. It is plausible that the (flat) ACDM cosmological 
model is breaking down. What remains is to confirm this diagnosis. 

To that end, a simple insight comes directly to us from the 
Friedmann equations; Ho is by construction an integration constant. 
In other words, it is theoretically a constant within the Friedmann- 
Lemaitre—Robertson—Walker (FLRW) framework on which ACDM 
is established by assuming spatial flatness. Nevertheless, it is obser- 
vationally a constant only when one has the correct cosmological 
model Krishnan et al. (2021), Krishnan & Mondol (2022). This 
statement is true not only for Ho, but also for other model parameters 
that are integration constants, such as matter density parameter 
(today) Qm := Pmo/3 He (Omo being the matter energy density today) 
in the ACDM model. Thus, if the ACDM model is breaking down, 
as all cosmological models must at some precision for the reason 
outlined earlier, one expects Hp, Qm, etc., to evolve with effective 
redshift. There are now numerous observations suggesting that Ho 
evolves, more accurately decreases with effective redshift Wong et al. 
(2020), Millon et al. (2020), Krishnan et al. (2020), Dainotti et al. 
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(2021), Dainotti et al. (2022), Ó Colgáin et al. (2022), Ó Colgáin et al. 
(2022), Malekjani et al. (2023) (see also Hu & Wang (2022), Jia, Hu 
& Wang (2023)). Likewise there are claims of Q,, increasing with 
effective redshift Ó Colgdin et al. (2022), Ó Colgdin et al. (2022), 
Malekjani et al. (2023) (see also Risaliti & Lusso (2019), Lusso et al. 
(2020), Yang, Banerjee & Ó Colgáin (2020), Khadka & Ratra (2020), 
Khadka & Ratra (2021), Khadka & Ratra (2022), Pourojaghi, Zabihi 
& Malekjani (2022), Pastén & Cárdenas (2023)). Furthermore, it 
has been noted in Ó Colgáin et al. (2022), Ó Colgáin et al. (2022), 
Malekjani et al. (2023) that evolution of Ho and Qm within ACDM are 
anticorrelated with each other. Note that if Q, evolves with effective 
redshift, it is unlikely that Sg is a constant, because og conspiring to 
balance evolution in Qm represents a contrived or unnatural scenario. 

Here we focus on the Sg discrepancy Aghanim et al. (2020), 
Heymans et al. (2013), Joudaki et al. (2017), Troxel et al. (2018), 
Hikage et al. (2019), Asgari et al. (2021), Abbott et al. (2022). Taken 
at face value, this anomaly says that Sg as inferred by the Planck 
collaboration (high-redshift inference) has a larger value than galaxy 
surveys (low-redshift measurement). This can be independently 
verified with growth rate data, in particular measurements of fog 
(z), which are independent of galaxy bias. Concretely, we will 
assume expressions that are valid for ACDM cosmology and fit 
fos (z) constraints in a single bin of increasing effective redshift. 
However, since growth rate data suffers from large fractional errors, 
we will leverage a working assumption in modern cosmology that 
Qn is tightly constrained to elucidate the trend. This means that the 
increases observed in Sg are driven by increases in og. Interestingly, 
this effect is also evident when one compares og constraints from 
number counts of galaxy clusters identified through the Sunyaev— 
Zeldovich effect at low redshifts with og constraints from Lya 
spectra at high redshifts Esposito et al. (2022). Moreover, CMB 
lensing, an observable most sensitive to redshift ranges z € [0.5, 5] 
and peaking in sensitivity at z ~ 2, recovers the Planck result Qu et al. 
(2023), Madhavacheril et al. (2023). This seemingly constrains any 
‘Sg tension’ to the late universe Marques et al. (2023).'! Admittedly, 
Qn may not be a constant Ó Colgáin et al. (2022), Ó Colgáin et al. 
(2022), but this is currently a fringe viewpoint. While Sg or øg tension 
in growth rate data is well studied in the literature Macaulay, Wehus & 
Eriksen (2013), Battye, Charnock & Moss (2015), Nesseris, Pantazis 
& Perivolaropoulos (2017), Kazantzidis & Perivolaropoulos (2018), 
Skara & Perivolaropoulos (2020), Quelle & Maroto (2020), Li et al. 
(2021), Benisty (2021), Nunes & Vagnozzi (2021), our goal in this 
letter is to explore redshift evolution of Sg within ACDM in the late 
universe under standard assumptions. 


2 WARM UP 


To get oriented, we impose a Gaussian prior, Qm = 0.3111 + 0.0056, 
which arises from combining Cosmic Microwave Background 
(CMB), galaxy, quasar, and Ly æ Baryon Acoustic Oscillation (BAO) 
constraints Aghanim et al. (2020), Alam et al. (2017), Hou et al. 
(2020), Neveux et al. (2020). In line with standard practice, we 
assume there is no discrepancy between CMB and BAO on 2,, 
inferences, that is, that the key success of the ACDM model is not 
undermined. The CMB+BAO prior is needed to compensate for 
the relatively low quality of the growth rate data, which will be 
further reduced by the removal of low-redshift data. Compared to 


‘Intriguingly, higher redshift observables may prefer lower values of og 
relative to Planck Miyatake et al. (2022), Alonso et al. (2023). Qm is poorly 
constrained at high redshifts. 
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Figure 1. fog (z) constraints in red with best fit values of (Qm, Sg) both with 
(blue) and without (green) a Gaussian prior Qm = 0.3111 + 0.0056. 


Planck Aghanim et al. (2020), our constraint is marginally more 
constraining, but remains representative of a Planck prior. In the 
absence of Gaussian priors, we employ uninformative uniform priors 
Qn ~ UO, 1) and og ~ U(0, 1.5) throughout. We combine this prior 
with 20 measurements of fos (z) from peculiar velocity and redshift- 
space distortion (RSD) data Said et al. (2020), Beutler et al. (2012), 
Huterer et al. (2017), Boruah, Hudson & Lavaux (2020), Turner, 
Blake & Ruggeri (2022), Blake et al. (2011), Blake et al. (2013), 
Howlett et al. (2015), Okumura et al. (2016), Pezzotta et al. (2017), 
Alam et al. (2021), as compiled recently in Nguyen, Huterer & Wen 
(2023). Modulo the removal of low-redshift data, we follow the 
methodology of Nguyen, Huterer & Wen (2023). We illustrate the 
constraints in Fig. 1 in red. We will be interested in the combination 


Ss [= Ogy Qn /0.3, a) 


where the constant og is the amplitude of matter fluctuations 
in spheres of 8 h~! Mpc with h := Ho/100kms~'Mpc™!. Follow- 
ing Wang & Steinhardt (1998), we introduce the matter density 
parameter, 

Qm (z) Qm(1 + 2)? 


SEE TOPI = Gn eo, 2) 


thereby allowing us to obtain 


6 z 2ú’) , 
fos(z) = og Q(z) exp @ a) : (3) 
0 


142’ 


As explained in Wang & Steinhardt (1998), these expressions are 
a valid approximation for ACDM. Concretely, we checked that the 
fractional error between the approximation and the exact expression 
for fog (z) based on the hypergeometric function F; (see Nesseris, 
Pantazis & Perivolaropoulos (2017)) is greatest at 1 percent at 
z = 0. In the range 0.35 S z <S 2, the fractional error is less 
than 0.2 per cent. This uncertainty is negligible compared to the 
observational uncertainties in Fig. 1. Working with an approximation 
may seem obsolete, but it allows us to shed light on evolution that 
is hidden in the analysis of Nguyen, Huterer & Wen (2023),” where 
the same approximation and data are employed. For all redshifts z, 
the function Q (z) is bounded Q,, < Q (z) < 1. In particular, Q 
— l as z — œ, so that Q,, dependence drops out as an overall 
factor in equation (3) at high redshifts, but some knowledge of Qm 
is retained through the integral. Thus, one expects Q,, to be poorly 


2In Nguyen, Huterer & Wen (2023) it is assumed that there is no evolution in 
cosmological parameters across the fog (z) constraints. 
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Figure 2. Same as Fig. | with Gaussian prior Qm = 0.3111 + 0.0056 and 
data with z < zmin removed. Sg increases with effective redshift. 


constrained at higher redshifts. As we show in the appendix, one 
risks encountering a degeneracy between the fitting parameters (Q, 
os) that is difficult to cleanly break with exclusively high-redshift 
binned data. Mathematically, one is guaranteed to run into a problem 
constraining Qm in high-redshift bins, so it is prudent to impose 
an Qn prior. Further, we will relax the prior and comment on the 
changes. 

The green curve in Fig. 1 denotes the best fit of the ACDM model 
to the 20 data points, whereas the blue curve documents the effect 
of introducing the Gaussian prior on Qm. It is worth noting that this 
increases the Sg inference and visibly worsens the fit to the fog (z) 
constraints at z © 0 and z © 1.5, thereby underscoring the tension 
between the Q, prior and the lower value of Q,, preferred by fo (z) 
data. The blue curve worsens the fit to fog (z) data by Ax? ~ 4.9. 
The fog (z) constraints are weak while the Gaussian prior is strong. 
Therefore, the prior effectively fixes Qm, so that the only parameter 
being fitted is og. As is clear from equation (3), one is simply fitting 
the scale of the fog (z) function, whereas the functional form is fixed. 
In Fig. 2, we see the effect of removing fog (z) constraints below 
Z = Zmin. We see that with Q,, constrained through the prior, og 
increases leading to larger values of Sg. As Sg increases, the curves 
visibly provide a better fit to high-redshift fo (z) constraints. This 
increase in Sg is driven by the two high-redshift data points and is 
expected, since it is visible in the raw data, that is, no data analysis 
is required. Note, the eBOSS data point at z ~ 1.5 is re-analysed in 
Brieden, Gil-Marin & Verde (2022) with different methodology and 
the central point evolves little. There is nothing to suggest this data 
point is not robust. 

We next split the fog (z) data at z = 0.7 into a low-redshift 
sample of 14 data points and a high-redshift sample of 6 data points. 
This choice is arbitrary, but it is guided by Fig. 1 and the intuition 
that evolution in Sg is expected once we see evolution in Hy) Wong 
et al. (2020), Millon et al. (2020), Krishnan et al. (2020), Dainotti 
et al. (2021), Dainotti et al. (2022), Ó Colgáin et al. (2022), Ó 
Colgáin et al. (2022), Malekjani et al. (2023) and Q,, Ó Colgáin 
et al. (2022), Ó Colgáin et al. (2022), Malekjani et al. (2023) 
elsewhere. We impose the Gaussian prior on Qm to each sample 
and perform a Markov Chain Monte Carlo (MCMC) analysis. The 
result is shown in Fig. 3, where in line with expectations from Fig. 2, 
the contours separate in the og or Sg direction. The contours are 
Gaussian, allowing us to directly compare Sg = 0.753705 (z < 0.7) 
to Sg = 0.839793) (z > 0.7) and conclude that the discrepancy in 
Sg is 1.60. Furthermore, the former value is discrepant with the 
Planck value Sg = 0.832 + 0.013 Aghanim et al. (2020) at 3.20. 
Throughout, we quote lo confidence intervals. Note that the 1.60 
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Figure 3. Marginalization of parameters (Qm, og, Sg) with Gaussian prior 
Qm = 0.3111 + 0.0056 and compilation of fog (z) data split at z = 0.7. Sg is 
reconstructed from (Qm, 03) MCMC chains. The discrepancy in the Sg plane 
between low- and high-redshift subsamples is 1.60. 


EE z2>07 

(l =m z<07 
+ + t +— 

i) 
4b 4 
12 4 
© 

8 10l ] 
08 p 4 
osr 4 È 4 È 1 
15 n gi 

02 04 06 08 06 08 10 12 14 0.5 1.0 15 

Qn Og Ss 


Figure 4. Same as Fig. 3 but without a Qn prior. 


shift occurs within the 20 data point sample, which simply assumes 
that BAO and CMB can consistently constrain Qm. Observe also that 
the effect of the shift is to remove a ~3o tension with Planck in the 
Sg parameter. 

In Fig. 4, we investigate the effect of removing the Qm prior. The 
low-redshift sample leads to (Qm, og) values consistent within lo 
with Planck, (Qm, og) = (0.2577}-02%, 0.803+0:088), however Sg = 
0.74470 :)99 is discrepant at 3.40. We emphasize that we sample 
(Qn, Og), but reconstruct Sg as a derived parameter from the 
MCMC chains. The higher redshift contours are shifted away from 
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Figure 5. The CMB+BAO prior Qm = 0.3111 + 0.0056 intersects the 
confidence intervals at larger values of Q,, resulting in a larger value of Sg 
in the high-redshift (z > 1.1) sample. The prior is consistent with low- and 
high-redshift constraints within 20, so constraints can be safely combined. 


their Planck values by more than lo throughout, (Qm, 08, Ss) = 
(0.12478 157, 0.9130 gay, 0.5859 191)» but remain consistent within 
2o. It is clear that lower values of Sg are driven by lower than 
expected values of Q,, in both low- and high-redshift subsamples. 
Note that caution is required when interpreting these shifts in the 
absence of an informative prior, as the high-redshift sample is small 
with only 6 data points. In Fig. 5, we demonstrate the consistency 
between the Qm prior and the confidence intervals from low- and 
high-redshift samples in the (Qm, Sg)-plane. The intersection of the 
prior lo confidence interval explains the preference for the higher 
Ss value from the high-redshift subsample. Since all constraints are 
within 20, no objection to combining the constraints is foreseen. 
While removing the Q, prior makes a considerable difference, it is 
nevertheless true that we see less tension with Planck in the high- 
redshift sample following shifts in cosmological parameters. 

As an aside, with an Q,, prior, we note that the ACDM model 
is still preferred over any theoretically poorly motivated model with 
a jump in og at z = 0.7. The reduction in the minimum of the x’, 
A Xin = —2.59 is not enough to overcome the penalty of introducing 
two additional parameters, an additional og and a redshift for the 
jump Zjump, even in the Akaike Information Criterion (AIC). Our 
observation here is essentially a microcosm of the ACDM tensions 
debate; despite seeing symptoms of a problem, here evolution in 
Sg when fog (z) constraints are combined with CMB+BAO, it is a 
separate matter to produce a model that outperforms ACDM. 

Given the relatively small size of the fog (z) data set considered in 
this work, in the appendix we work with a larger data set of historical 
data. The shortcomings of this data set are explicitly discussed, but 
we arrive at similar conclusions. 


3 CONCLUSIONS 


Either from the Friedmann equations directly, or the continuity equa- 
tion and the assumption of pressure-less matter, both Ho and Q,,, arise 
mathematically as integration constants. Thus, consistency demands 
that both Hp and Q,, are observationally constants. An important 
corollary is that CMB, BAO, Type Ia supernovae, etc., must provide 
consistent Qm constraints in ACDM cosmology. Similarly, the fitting 
parameter og must be a constant, as it is related to the normalization 
of the matter power spectrum. Here, working within the assumption 
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that CMB and BAO consistently constrain Q,,, we have shown that 
the Sg := og./Q,,/0.3 parameter increases with effective redshift 
in the late universe, z S 2. This evolution is compelling because 
it corroborates the tendency of cosmic shear to give lower values 
of Sg than Planck-ACDM Aghanim et al. (2020), Heymans et al. 
(2013), Joudaki et al. (2017), Troxel et al. (2018), Hikage et al. 
(2019), Asgari et al. (2021), Abbott et al. (2022). Moreover, it is 
not an isolated observation; we note that (i) og constraints from 
low-redshift clusters are biased lower than high-redshift Ly œ spectra 
Esposito et al. (2022) and (ii) weak-lensing Sg constraints, sensitive 
to lower redshifts, are biased lower than CMB lensing results Qu et al. 
(2023), Madhavacheril et al. (2023), Marques et al. (2023) that are 
sensitive to higher redshifts. In Esposito et al. (2022), Qu et al. (2023), 
Madhavacheril et al. (2023), Marques et al. (2023), systematics are 
a greater concern as one is comparing different observables, but 
here we are working with common fog (z) constraints throughout, 
so systematics should be under greater control. The data largely 
comprises RSD but peculiar velocity constraints may be present at 
z ~ 0, for example Boruah, Hudson & Lavaux (2020), Said et al. 
(2020). However, these differences aside, we agree on the increasing 
og trend, a trend that is evident from the raw data. 
It is instructive to recall the assumptions being made: 


(i) Equation (3) is a valid approximation for ACDM behaviour. 

Gi) Qm in ACDM is approximately 0.3. 

(iii) Data sets and/or priors from independent data sets can be 
combined. 

(iv) The data are correct. 


If this set of assumptions is correct, then the increasing o/S¢ 
trend follows. Assumption (i) is standard Wang & Steinhardt (1998) 
and the theoretical uncertainty introduced is less than 1 per cent, 
more accurately 0.2 per cent in the redshift range 0.35 S z S 2. 
Assumption (ii) is a widely recognized observation. Assumption (iii) 
is standard practice. Assumption (iv) is the weakest, but the removal 
of data points requires justification. In other words, removing 
data points with the sole motivation of recovering a null result is 
unscientific. 

Concretely, we presented analysis in a sample of 20 and 66 fog 
(z) (see appendix for the latter) constraints respectively. In both 
samples we find the same increasing Sg trend, which is essentially 
a og trend, because Q, is subject to a strong prior. Since there 
is considerable survey overlap in Appendix Fig. A3, the tension 
between low- and high-redshift is admittedly overestimated at 2.80 
( Appendix Fig. A3). That being said, if the 20 data points in Fig. 1 
are independent (Nguyen, Huterer & Wen (2023) assumes they are), 
then there is a 1.60 shift in Sg that warrants further study, especially 
since it seems to ameliorate tension with Planck. One could identify 
other redshift ranges where the discrepancy between low- and high- 
redshift inferences is less, but such an exercise is meaningless. To 
stress test any sample against evolution, one needs to focus on splits 
that exacerbate the feature. Moreover, Appendix Fig. A2 clearly 
demonstrates that evolution is present throughout the sample. Despite 
the presence of evolution, we find that the ACDM model is still 
preferred over any theoretically ad hoc model with a jump in Sg. 
Finally, when removing the Q, prior, despite noticeable differences 
in cosmological parameters, we find that high-redshift data are more 
consistent with Planck. 

Forthcoming releases from the Dark Energy Spectroscopic In- 
strument (DESI) DESI Collaboration et al. (2016) will be in a 
position to confirm or refute our assumption (iv). Thus, our result 
is preliminary. Our findings also provide a targeted prediction for 
tomographic studies of cosmic shear Hikage et al. (2019), Asgari 
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et al. (2021), Abbott et al. (2022). These surveys, that is, HSC, 
KiDS, and DES, currently employ redshift bins with varying effective 
redshifts across the surveys and shifts in cosmological parameters 
up to lo are evident. Given our results here, it is plausible that 
Sg redshift evolution will be detected as cosmic shear data quality 
improves through a better understanding of systematics. If redshift 
evolution is not detected, and assuming Sg tension is physical, then 
there must be evolution with scale Marques et al. (2023). Finally, 
our results caution theoretically that physics that alters the radius of 
the sound horizon Knox & Millea (2020) cannot account for these 
expected hallmarks of model breakdown. 
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APPENDIX A: LARGER ARCHIVAL DATA SET 


Given that the increase in Sg in the text is driven exclusively by two 
high-redshift data points, it is prudent to work with a larger data set. 
The role of the second data set is to provide a sanity check, since the 
larger data set comes with caveats that we discuss further. To that 
end, we focus on table II of Kazantzidis & Perivolaropoulos (2018), 
where one finds 63 historical measurements of fog (z) Beutler et al. 
(2012), Blake et al. (2013), Howlett et al. (2015), Okumura et al. 
(2016), Huterer et al. (2017), Song & Percival (2009), Blake et al. 
(2012), Davis et al. (2011), Hudson & Turnbull (2013), Turnbull 
et al. (2012), Samushia, Percival & Raccanelli (2012), Tojeiro et al. 
(2012), de la Torre et al. (2013), Chuang & Wang (2013), Sanchez 
et al. (2014), Wang et al. (2018), Feix, Nusser & Branchini (2015), 
Chuang et al. (2016), Alam et al. (2017), Beutler et al. (2017), Wilson 
(2017), Gil-Marin et al. (2017), Hawken et al. (2017), de la Torre 
et al. (2017), Pezzotta et al. (2017), Feix, Branchini & Nusser (2017), 
Howlett et al. (2017), Mohammad et al. (2018), Shi et al. (2018), 
Gil-Marin et al. (2018), Hou et al. (2018), Zhao et al. (2019). We 
add more recent data from Bautista et al. (2020), de Mattia et al. 
(2021), Neveux et al. (2020) to increase the sample to 66 data points. 
We observe that Refs. Beutler et al. (2012), Blake et al. (2013), 
Howlett et al. (2015), Okumura et al. (2016), Huterer et al. (2017) 
are common to both data sets in Fig. 1 and Fig. Al. The other 
data points are the work of independent groups, but as cautioned 
in Kazantzidis & Perivolaropoulos (2018), most analyses of growth 
rate data, e.g. Nguyen, Huterer & Wen (2023), tend to pick 20 odd 
data points to counteract any overcounting of data. However, this 
choice of independent data points is largely subjective. Moreover, 


Sg increases with redshift in ACDM 125 


— (2Zmin, Sg) = (0, 0.741) 
— (Zmin. Ss) = (0.3, 0.751) 
— (Zmin, Sg) = (0.7, 0.781) 
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Figure Al. fog (z) constraints alongside best fit values of Sg from fitting in 
equation (3) in the range z > Zmin. Throughout we employ the Gaussian prior 
Qm = 0.3111 + 0.0056. 
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Figure A2. Removing fog (z) constraints with z < Zmin while fitting 
expressions for the ACDM model in equation (3) subject to a Gaussian 
prior Qm = 0.3111 + 0.0056. An increasing Sg trend with effective redshift is 
evident and the error bars correspond to 1ø confidence intervals. Constraints 
from KiDS, DES, and Planck are illustrated for comparison. 


different studies make use of different fiducial cosmologies. To 
address this latter concern, we follow equation (2.6) of Nesseris, 
Pantazis & Perivolaropoulos (2017) (see also Macaulay, Wehus 
& Eriksen (2013)) and correct for this choice of fiducial model 
by performing a rescaling by the combination H(z)D,(z), where 
A(z) and D(z) denote the Hubble parameter and angular diameter 
distance at a given z. This means that the data set is consistent with 
the Q,, prior by construction. Note, following Nguyen, Huterer & 
Wen (2023), we did not perform this correction in the body of the 
letter. In the fog (z) literature, one finds groups who correct and do 
not correct for the fiducial cosmology. 

In Fig. A2 we show the effect of removing fog (z) constraints 
below z = Zmin. Where relevant, for example WiggleZ Blake et al. 
(2012), we crop covariance matrices accordingly as we remove data 
points. Our strong Gaussian prior on Q, ensures that our a posteriori 
distributions from MCMC analysis remain Gaussian. As a result, the 
lo errors in Fig. A2 are representative. The reader will note that, 
in line with expectations, the size of the error bar increases as we 
remove data. Furthermore, the error bars in Fig. A2 are correlated, 
because all plotted constraints larger than a given Zmin share data 
points. 

We can address criticism regarding the correlated constraints in 
Fig. A2 by simply splitting the sample into independent subsamples. 
Our results in Fig. A2 suggest that any split at z ~ 1 may maximize 
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Figure A3. Marginalization of parameters (Qm, og, Ss) with Gaussian prior 
Qm = 0.3111 + 0.0056 and compilation of fog (z) data split at z = 1.1. Sg is 
reconstructed from (Qm, 03) MCMC chains. The discrepancy in the Sg plane 
between low- and high-redshift subsamples is 2.80. 
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Figure A4. Same as Fig. A3 but without the CMB+BAO Qn prior. The 2D 
posterior in the (Qm, 7g)-plane points to a degeneracy that is poorly broken 
by higher redshift constraints since the lo confidence intervals saturate the 
bounds. 


the discrepancy between low-redshift and high-redshift subsamples. 
Once again, the trend is directly visible in the data in Fig. A1, since 7 
from 8 data points above z = 1 prefer larger values of S/o than the 
full sample (Zin = 0). The choice of split is somewhat arbitrary, as it 
depends on the composition or makeup of the overall sample, notably 
redshift distribution and quality of the data. For this reason, the split 
here is conducted at z = 1.1, whereas it was z = 0.7 in the body 


of the letter. Note that our sample here now has an extended range, 
so we can explore higher redshift splits of the sample. Concretely, 


we find that a split at z = 1.1 leads to Sg = 0.732 + 0.014 at lower 
redshifts and Sg = 0.908 + 0.062 at higher redshifts, with a tension 
at 2.80 as illustrated in Fig. A3. Once again, it is worth noting that 
Sg from data below z = 1.1 is in tension with Planck at ~5o, while 
the data above z = 1.1 is consistent with Planck at ~lo. 

Removing the Qm prior, the corresponding results 
are (Qm, Sg) = (0.229995! 0.823497) (z < 1.1) and 
(m, Sg) = (0.43579 397, 0.8947) 48) (Z = 1.1). Evidently, all 
inferences are now consistent with Planck within 20, however as 
can be seen from Fig. A4, there is degeneracy in the (Qm, o g)-plane 
that the data fails to break. Thus, care needs to be taken with these 
inferences as they will be prior dependent. Noting that Q, is lower 
than expected in the lower redshift subsample, it is interesting to 
studying the full sample without the Qm prior. Doing, so we find 
(Qm, og, Sg) = (0.19179-935, 1.11849 183, 0.890+0:066), Evidently 
there is some discrepancy with Planck, whereby Qm, og, and Sg are 
respectively 3.50, 1.90, and lo away. In both data sets considered 
in this work, there is evidently tension with Planck in one of the 
parameters. 

Bearing in mind that not all the data points are independent, and as 
a result the errors are underestimated, the main take away from our 
analysis here is that the Sg trend with Planck+-BAO Q, prior is robust 
to changes in the data set. A secondary take away is confirmation 
that one expects to run into trouble in higher redshift bins without an 
Qm prior in line with equations (2) and (3). 
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